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Abstract — Pulse position modulation (PPM) is a popular signal 
modulation technique which creates M-ary data by means of 
the position of a pulse within a time interval. While PPM 
and its variations have great advantages in many contexts, this 
type of modulation is vulnerable to loss of synchronization, 
potentially causing a severe error floor or throughput penalty 
even when little or no noise is assumed. Another disadvantage is 
that this type of modulation typically offers no error correction 
mechanism on its own, making them sensitive to intersymbol 
interference and environmental noise. In this paper we propose a 
coding theoretic variation of PPM that allows for efficient symbol 
and frame synchronization as well as strong error correction. The 
proposed scheme can be divided into a synchronization layer 
and a modulation layer. This makes our technique compatible 
with major existing techniques such as standard PPM, multipluse 
PPM, and expurgated PPM as well in that the scheme can be 
realized by adding a simple synchronization layer to one of 
these standard techniques. We also develop a generalization of 
expurgated PPM suited for the modulation layer of the proposed 
self -synchronizing modulation scheme. This generalized PPM can 
also be used as stand-alone error-correcting PPM with a larger 
number of available symbols. 

Index Terms — Pulse position modulation, PPM, synchroniza- 
tion, error correction, self-synchronizing code, comma-free code, 
combinatorial design, optical orthogonal code. 



I. Introduction 

THIS work addresses the problem of symbol and frame 
synchronization during signal modulation by means of 
the position of each pulse in the time domain when the 
channel is assumed to be noisy. The approach taken here is to 
rethink the problem in a coding theoretic framework and give 
a theoretical and general foundation. 

Among various forms of information, perhaps the simplest 
of them is binary information. In many communications sce- 
narios, however, it is beneficial to transmit data in M-ary 
format with M > 2. For this reason, there have been proposed 
various modulation techniques that support not only binary 
format but also M-ary format with large M. 

Pulse position modulation, or PPM for short, is one of the 
more popular modulation techniques [1|. In this modulation 
technique, each symbol occupies a time interval of equal 
length. A symbol interval is divided into Q time slots of 
equal length, where exactly one pulse is transmitted at one 
of the Q time slots. Which symbol each interval represents 
is determined by at which time slot the unique pulse is 
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transmitted. Because there are Q choices of pulse positions for 
each symbol interval, PPM offers M = Q distinct symbols. 

While PPM is one of the most fundamental forms of signal 
modulation in use today, there are some inherent drawbacks. 
For instance, in a communications system with a severe peak 
power constraint, PPM becomes inefficient as the number 
of symbols increases because the energy per symbol drops 
accordingly. 

Multipulse PPM is a generalization of PPM to mitigate this 
problem, where K pulses are sent during each symbol interval 
by using K out of Q time slots so that M = (^) symbols 

can be represented J2]. Let Vk — yVi € | wt(i>i) = K j 
be the set of all Q-dimensional binary vectors Vi of weight 
K. We let Is represent time slots at which single pulses are 
transmitted and Os those at which no pulse is sent. From 
the viewpoint of coding theory, the symbols of PPM can 
be seen as the binary constant-weight code V\ of length Q, 
weight one, and minimum distance two with Q codewords, 
whereas the symbols of multipluse PPM can be regarded as 
the binary constant-weight code Vk of length Q, weight K, 
and minimum distance two with (^) codewords. 

The non-coherent nature of PPM and multipulse PPM 
makes them attractive to communications systems in which 
coherent detection is expensive or impossible, such as typical 
optical communications [3|. However, as is evident from the 
fact that the codes V\ and Vk are both of minimum distance 
two, one of the major disadvantages of these modulation 
techniques is that they are inherently vulnerable to intersymbol 
interference and natural environmental noise. Hence, PPM and 
multipulse PPM require strong error correcting schemes at a 
higher level when errors are of concern. 

Expurgated PPM is a recently proposed modulation tech- 
nique that generalizes PPM in a way error correction can 
be provided at the modulation stage while offering the same 
number M = Q of symbols [4|. If we see it from the coding 
theoretic point of view, the key idea of expurgated PPM can 
be understood as using special combinatorial designs called 
symmetric designs to define constant-weight codes with large 
minimum distances which allows for simple implementation. 

The above three modulation techniques still share the 
other major disadvantage, namely the susceptibility to loss 
of synchronization. As in earlier research on synchronization 
for this type of modulation ll5l- lfT0l . we assume that slot 
synchronization is always provided so that the magnitude of 
misalignment can by expressed by a multiple of the length 
of a time slot. Under this assumption, erroneous symbol 
synchronization means that the window of the receiver is 
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aligned to the consecutive Q time slots that consist of the last 
i time slots of one symbol interval and the first Q — i time 
slots of the following one for some positive integer i < Q — 1. 
Erroneous frame synchronization is understood the same way 
by regarding a set of consecutive / symbol intervals as a frame 
of length / with fQ time slots. 

In PPM, even if the channel is completely noiseless, er- 
roneous symbol synchronization can not be detected if there 
happens to be exactly one pulse within the misaligned win- 
dow. Similarly, the receiver can not detect erroneous symbol 
synchronization under multipulse PPM if there are exactly K 
pulses within the misaligned window. In the case of expurgated 
PPM, erroneous symbol synchronization will go unnoticed if 
the Q slots in the misaligned window form a valid codeword 
or invalid one within the decodable range. Hence, without 
some sort of synchronization mechanism, these modulation 
techniques may exhibit a severe error floor. 

The known method for alleviating this synchronization 
problem in the literature is to periodically insert a synchroniza- 
tion marker that consists of sQ time slots Tj, < i < sQ — 1, 
for some positive integer s in which for any nonnegative 
integer j < s — 1 the Q consecutive time slots T i+ jQ, < i < 
Q — 1, form a valid codeword (see, for example, [5], [6|, [10|). 
In other words, a certain pattern of consecutive s symbols is 
periodically inserted to signal the boundaries. Ideally, the off- 
peak autocorrelations of the synchronization marker should 
be as small as possible to suppress the probability that the 
receiver misses the marker or is deceived by a false one 
at an unintended position. For PPM and multipulse PPM, it 
is possible to find a synchronization marker that causes no 
ambiguity as long as the channel is noiseless. However, an 
ambiguous synchronization marker can not be short. If the 
channel is noisy, synchronization with this approach becomes 
increasingly difficult and complicated as well. 

The primary purpose of this paper is to provide a unified 
solution to both error tolerance and synchronization by using 
coding theory. We propose a modulation scheme in which reli- 
able symbol and frame synchronization can be achieved under 
the presence of noise through the same type of correlation 
receiver used for expurgated PPM. The required overhead in 
terms of the size of a synchronization marker is significantly 
smaller than the known method. The signal modulation can be 
performed separately, so that the known modulation techniques 
of PPM kind can be exploited straightforwardly. We also 
develop a generalized format of expurgated PPM which has a 
larger number of symbols and increased minimum distance 
than standard PPM and expurgated PPM. This generalized 
PPM can also be used for the proposed self-synchronizing 
modulation scheme in place of other modulation techniques 
of PPM kind to offer higher error tolerance and/or better 
throughput. 

II. The scheme 

We describe our scheme as a special class of constant- 
weight codes by exploiting the coding theoretic view intro- 
duced in the previous section. For the sake of generality, a sim- 
ple setting is assumed where the probability that the receiver 



fails to correctly decode symbols decreases monotonously as 
the minimum distance increases. Thus, for the most part, our 
framework will focus on the minimum distance and generally 
aim for the largest possible codewords for given parameters. 

This section is divided into two subsections. Subsection lll-AI 
is devoted to our self-synchronization mechanism that allows 
for separate implementation of a modulation layer based on 
pulse positions in the time domain. Subsection IH-BI gives a 
generalized version of expurgated PPM that can be used both 
as a stand-alone modulation scheme with error correction and 
as the modulation part of our scheme. 

A. Synchronization layer 

As is pointed out in the previous section, PPM and its 
variations can be interpreted in terms of constant-weight codes. 
In this coding theoretic framework, the number of codewords 
corresponds to the number of available symbols while the 
minimum distance corresponds to the error tolerance capabil- 
ity. The constraint that the weight of a codeword is constant 
ensures equal energy across symbols. Hence, a good signal 
modulation technique based on the pulse positions within 
symbol intervals of fixed length corresponds to a constant- 
weight code of large minimum distance with many codewords 
that can be decoded in a certain simple manner. 

Since this coding theoretic view offers a clear picture of how 
the sensitivity to noise may be alleviated, it would be natural to 
ask if it is also possible to exploit the framework to overcome 
the other major weakness, namely the vulnerability to loss 
of synchronization, while preserving other major features of 
existing standard techniques. This subsection answers this 
question in the affirmative by developing a theory of symbol 
and frame synchronization for signal modulation of PPM 
type. In what follows, we use binary codes to represent 
symbols expressed by pulse positions, where Os in a codeword 
correspond to time slots with no pulse while Is represent those 
at which single pulses are sent. 

A self-synchronizing code C C F!, 1 is a binary block code of 
length n where the symbol string formed by an overlapped 
portion of any two concatenated codewords is not a valid 
codeword. In the coding theory literature, self-synchronizing 
codes are also called comma-free codes. The property that no 
codeword appears as a substring of two adjacent codewords 
allows for block synchronization without any external help as 
long as synchronization is provided at the bit level. 

The key idea of our approach is to use a special self- 
synchronizing code that allows for synchronization by observ- 
ing only part of the window and frees up the rest for modula- 
tion. To obtain constant-weight self-synchronizing codes with 
desirable properties for our purpose, we employ combinatorial 
design theory. 

Take a sequence of codewords of a binary block code 
of length n. A splice of length n between codeword 
) and the following codeword y = 
(yo, 2/i j • ■ • , Vn—x) m me codeword sequence is a concatenated 
binary sequence (x n -i, . . . , x n -x, y , . • . y n -i-i) composed 
of the last i bits of x and the first n — i bits of y for some 
positive integer i < n — 1. 
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A binary block code of length n is said to be of comma- 
free index p if the Hamming distance between any codeword 
z and any splice of length n between any two codewords 
x, y is at least p. It is straightforward to see that with a hard- 
decision algorithm a self-synchronizing code of comma-free 
index p assures block synchronization under the presence of 
up to L^t^J bit flips (or errors) in the received message of 
length n. In the AWGN channel, for example, one may use a 
correlation receiver for soft-decision synchronization to take 
advantage of the Hamming distance between a valid codeword 
and a splice. 

A difference system of sets (DSS) of index p over Z n is a 
family of disjoint subsets Di of Z n such that the multi-set 

{a-b (mod n) \ a £ Dt,b £ Dj,i ^ j} (1) 

contains every d £ Z n \ {0} at least p times. The difference 
between two elements from different subsets of Z„ is called 
an outer difference. A DSS is perfect if the multi-set defined 
in (Q]i contains every d £ Z n \ {0} exactly p times. A DSS 
is regular if all subsets Di are of the same size. For instance, 
the set {{1, 2}, {3, 5}} over Zg forms a regular DSS of index 
one because every nonzero outer difference appears at least 
once as follows: 

1-3 = 6 (mod 8), 3-1 = 2 (mod 8), 

1- 5 = 4 (mod 8), 5-1 = 4 (mod 8), 

2- 3 = 7 (mod 8), 3-2=1 (mod 8), 
2-5 = 5 (mod 8), 5-2 = 3 (mod 8). 

This DSS is not perfect because 4 appears twice as an outer 
difference while each of the other elements of Zg \ {0} occurs 
exactly once. 

The original motivation of the study of DSSs was to realize 
self-synchronizing codes as cosets of linear codes in order 
to achieve low encoding and decoding complexity QT), |12|. 
However, DSSs appear to have far greater potential and can 
be exploited to provide self-synchronizing codes with various 
desired properties. For our purpose, we would like binary 
constant-weight self-synchronizing codes of sufficiently large 
comma-free index with a cartain additional property. We use 
DSSs with exactly two sets to obtain desirable codes. 

Theorem 2.1: If there exist a DSS {Do, D\] of index p over 
Z„ and a binary constant- weight code C of length n — \Dq \ — 
\D\\ and weight K, then there exists a binary constant-weight 
self-synchronizing code of length n, comma-free index p, and 
weight K + \D\\ with \C\ codewords. 

Proof: Let {Do, Di} be a DSS of index p over Z„ and 
C a binary constant-weight code of length n — |Z?o| — \D\\ 
and weight K. For every codeword c £ C, construct the n- 
dimensional vector d c = (do, . . . , d„_i), where = for all 
i € D , d t — 1 for all i £ D-y, and the (n - \D \ - |Di|)- 
dimensional vector (di), i £ Do U D\, forms c. Each of the 
resulting \C\ elements of F2 is of weight K + \D\\. It suffices 
to show that the set T> = {d c | c £ C} of these vectors forms 
a code of comma-free index p. Take a pair of not necessarily 
distinct codewords from T> and form a splice s of length n 
by concatenating the first s bits of one codeword and the last 
n — s bits of the other for a positive integer s < n — 1. Take 



another not necessarily distinct codeword d from T>. Because 
{Do, Di} forms a DSS of index p, there are at least p ordered 
pairs (a, b) such that a — b = s (mod n), where a and b 
belong to different sets of the DSS. Thus, there are at least 
p discrepancies between s and d within the coordinates i £ 
flo U Di. The proof is complete. ■ 
To make the virtue of the above construction clearer, take 
the set {{1, 2, 3, 4, 5}, {0, 6, 11, 16, 21}}, which forms a per- 
fect regular DSS of index two over Z26- The two sets of 
cardinality five specify the positions of 0s and Is respectively 
as synchronization markers while the remaining 16 positions 
are freely available for signal modulation by a constant-weight 
code of length 16 such as multipulse PPM with Q = 16. If 
we use {1,2,3,4,5} for 0s and {0,6,11,16,21} for Is, by 
writing a bit used by the constant- weight code of length 16 as 
*, we have 26 bit sequence 

1000001****1****1****1****. 

Because each nonzero outer difference appears twice in the 
DSS, regardless of the content of each *, there are at 
least two discrepancies among the positions {1,2,3,4,5} U 
{0,6,11,16,21} between any pair of a valid codeword of 
the resulting self-synchronizing code and a splice. To give 
a smaller example, we can combine the DSS {{1, 2}, {3, 5}} 
over Zg and a constant-weight code of length four such as PPM 
with Q = 4 in the same way. In this case, the coordinates 0, 
4, 6, and 7 correspond to free bits. Note that in each of these 
examples the cardinalities |Do| and |Z?i| are the same. Hence, 
we may swap the roles of the sets to obtain a binary constant- 
weight self-synchronizing code of the same parameters with 
the same number of free bits. 

A self-synchronizing code constructed by the method given 
in Theorem 12.11 can be synchronized by only looking at 
periodic autocorrelations of the partial window specified by 
the corresponding DSS. Hence, the synchronization device 
on the receiver side only needs the corresponding signals as 
inputs for synchronization. If we use symbol intervals of PPM, 
multipulse PPM or expurgated PPM as the constant-weight 
code C of length Q = n — \Do\ — \D\\ in Theorem 12.11 
we achieve symbol synchronization that is securely checked 
for each symbol interval while allowing for modulation by 
pulse positions in the time domain through freely available 
bits. Packing / symbol intervals into the freely available part 
provides frame synchronization that is constantly checked for 
each frame consisting of / symbol intervals. If we use a DSS 
over Z„ for symbol synchronization and then employ another 
DSS {D'o, D[} overZ. n/ such that n'- \D' a \- \D[ \ is a multiple 
/ of n, then the scheme provides both symbol and frame 
synchronization. 

In the remainder of this subsection, we explore properties 
of DSSs and demonstrate that our synchronization method is 
significantly more efficient than the known technique that uses 
a sequence of symbols as a synchronization marker. 

Because a DSS {Do,D{\ over Z n leaves n — \Do\ — \Dy\ 
bits for the modulation layer, ceteris paribus, it is desirable 
for |Z?o| + l^il to be small. This parameter is called the 
redundancy of a DSS. In our context, redundancy is the 
parameter that denotes the number of time slots we sacrifice 



4 



PREPRINT TYPSET IN IEEE TRANSACTIONS STYLE 



for synchronization per symbol interval or per frame interval. 
We use r(n, p) to denote the smallest achievable redundancy 
for given order n and index p. A DSS is optimal if its 
redundancy is r(n,p). The following is a special case of the 
well-known Levenshtein bound: 

Theorem 2.2 (fiTTH): For any DSS of index p over Z„ with 
exactly two sets, it holds that 



r{n,p) > y/2p{r 



1) 



with equality if and only if the DSS is perfect and regular. 

If we are allowed to have a sufficiently strong signal for 
each pulse compared to noise, it may be enough to employ a 
DSS of index one or two. In this case, the following classical 
results give optimal DSSs: 

Theorem 2.3 { /[TTil . ^T3^): For any integer n > 2, the pair 
of sets 

Do = {ir + 1 | 1 < i < n} and 
Di = {i\l<i< t } 

form an optimal DSS of index one over Z„, where 



T 



1 



2r 



and Ti 



1 



Note that when the order of the ring Z„ is of the form n = 
2m 2 + 1 for some positive integer m, the redundancy of the 
optimal DSS given above achieves the Levenshtein bound with 
equality. 

Theorem 2.4 (fiTTH): For any integer n > 2, the pair of sets 

D = {t + 1} U {n - ito | < i < n - 2} and 
£>i = {i | 1 < i < t } 
form an optimal DSS of index two over Z n , where 
"n- 1" 



T() 



T() 



and n — \y/n — 1 1 



As in Theorem 12.31 the redundancies of the optimal DSSs in 
Theorem 12.41 achieve the lower bound in Theorem 12.21 with 
equality when n = m 2 + 1 for some positive integer m. 

It is notable that these optimal DSSs only require 0(n^) 
bits for synchronization. For instance, if we use an opti- 
mal DSS for symbol synchronization with PPM, because 
n = 0(Q), the number of time slots we sacrifice is 0(Qz) 
per symbol interval. Any method that inserts a sequence of 
valid symbol intervals must sacrifice sQ — ft(Q) time slots 
for some positive integer s, occupying a significantly larger 
number of time slots than our method. It is also worth noting 
that Theorems 12.31 and 12.41 explicitly give optimal examples 
for all nontrivial order n > 2. 

If the peak power is severely limited, we may need a DSS 
of larger index for secure synchronization. Such DSSs have 
been studied in various contexts (see, for example, lfl4l - 
[22 1 and references therein). To study the use of a DSS 
{-Do, Di} over Z„ for synchronization in the context of signal 
modulation by pulse positions, we use the redundancy rate 
R = \ D °\ + \ Dl \ of the DSS to measure its slot usage. For 
instance, if the redundancy rate is half, synchronization and 
modulation require the same amount of time resources. The 



least useful DSSs for our purpose are those of redundancy 
rate one because there would be no time slots for modulation. 
We aim for the smallest possible redundancy rate for given 
n and p. Note that if we define an analogous parameter for 
a synchronization method that inserts a sequence of valid 
codewords for symbol synchronization by taking the fraction 
between the numbers of time slots for synchronization and 
those for modulation per symbol, such a value can never be 
less than a half because one ought to insert at least one symbol 
per symbol interval. As we will see next, our method can break 
this fundamental limit even when the channel is assumed to 
be noisy. 

The lower bound on the achievable redundancy given in 
Theorem 12.21 suggests that the number of time slots required 
for synchronization may still be only 0(n^) for the case when 
the index is larger than two. In fact, there are known classes of 
optimal DSSs with exactly two sets that meet the Levenshtein 
bound with equality. In the remainder of this subsection, we 
list known optimal DSSs that are useful for our purpose as 
well as almost optimal DSSs that are equally of interest but 
are not found in the literature. 

Let p = eg + 1 be an odd prime for some positive integers 
e and g. The eth cyclotomic classes in ¥ p are defined as Cf = 
{a l+te | < t < g — 1}, where a is a primitive element of ¥ p 
and < i < e — 1. The cyclotomic numbers (i,j) e of order e 
are defined as (i,j) e = \ {Cf + 1) U Cf\. Note that the ring Z p 
may be regarded as the finite field F p by defining the natural 
division when p is prime. We use the following special case 
of the construction for DSSs given in [23 1: 

Theorem 2.5 ( R23\l ): Let p — 2eh + 1 be an odd prime, 
where e and h are positive integers. The set {Cg e ,Cg e } of 
two cyclotomic classes in ¥ p forms a regular DSS over Z p 
and index p, where 

p = min{(z, e) 2e + (i + e, e) 2e | < i < e - 1} . 



In particular, if 



(«, e) 2e + + e, e) 2e 



for every i, then the regular DSS is of index -, perfect, and 
hence optimal. 

The following are the two known classes of optimal DSSs 
with exactly two sets constructed in this manner: 

Theorem 2.6 ( H4\l ): For every m such that 16m 2 + 1 is 
an odd prime, the set {Cq , C|} of two cyclotomic classes 
in F 16m 2 +1 forms a perfect regular DSS of index 2m 2 and 
redundancy rate \ - 32m ^ +2 over Z 16m 2 +1 . 

Theorem 2.7 ( Hl4\l ): For every m such that 108m 2 + 1 is 
an odd prime, the set {Cq,C®} of two cyclotomic classes 
in F 10 8 TO 2 +1 forms a perfect regular DSS of index Qm 2 and 
redundancy rate | - 324 ^ +3 over Z 108m 2 +1 . 

Because the above optimal DSSs are both perfect and regu- 
lar at the same time, their redundancies meet the Levenshtein 
bound with equality. For instance, the DSS over Z 16m 2 +1 in 
Theorem 12.61 uses 8m 2 time slots for synchronization and 
leaves the remaining 8m 2 + 1 time slots for modulation. 
For comparison to the method that inserts a sequence of 
symbols, in the case of symbol synchronization, this means 



FUJIWARA: SELF-SYNCHRONIZING PULSE POSITION MODULATION WITH ERROR TOLERANCE 



5 



TABLE I 

Perfect regular DSSs from cyclotomic constructions for 
m < 10 



3 + 2x 



m 


n 


l-Dfl 


P 


— \ D a l + l-Di 

n 


Reference 


1 


17 


4 


2 


8 

17 


Theorem |2.6| 


1 


257 


64 


32 


128 
257 


Theorem 12.61 


5 


401 


100 


50 


200 
40.1 


Theorem 12.61 


6 


577 


144 


72 


2 s s 
577 


Theorem 12.61 


9 


1297 


324 


162 


048 
12' >7 


Theorem 12.61 


10 


1601 


400 


200 


Si in 

1601 


Theorem 12.61 


1 


109 


18 


6 


36 
10!) 


Theorem 12.71 


2 


433 


72 


24 


144 
455 


Theorem 12.71 


6 


3889 


648 


216 


1200 

3889 


Theorem |2.7| 



that the efficiency of the DSSs in Theorem 12.61 in terms of 
slot usage is almost the same as inserting only one symbol per 
symbol interval. Any method that inserts a symbol sequence 
would not be able to assure synchronization at this redundancy 
rate even if the channel is almost noiseless, whereas our 



method tolerates up to |_- 



i 2 — 1 bit flips for hard- 
decision decoding and an equivalent level of noise for soft- 
decision decoding. By the same token, Theorem 12.71 gives 
DSSs in which the number of time slots per symbol interval 
for synchronization is only about a half of the number of time 
slots for modulation. As their redundancy rate R~\ suggests, 
this level of efficiency in terms of the number of sacrificed 
time slots is fundamentally unachievable by any method that 
inserts valid codewords for synchronization. Table [I] lists the 
perfect regular DSSs obtained by Theorems 12.61 and 12.71 for 

TO < 10. 

If we allow the redundancy of a DSS to be slightly above 
the right-hand side of the lower bound in Theorem 12.21 we 
may extend the range of parameters covered by the same 
construction method while keeping redundancy very low. This 
idea was investigated in a more general setting by Mutoh in 
his unpublished manuscript [24|. Here we give two classes 
of DSSs by simply plugging cyclotomic numbers calculated 
in ll25ll into Theorem 12.51 For details of the calculations of 
cyclotomic numbers, we refer the reader to 11251 . ll26l . 

Theorem 2.8: Let n = 1 (mod 4) be a prime and binary 
quadratic form n = x + Ay 2 with x = 1 (mod 4) its 
decomposition. Then the set {CqjC^} of two cyclotomic 
classes in F n forms a regular DSS of index p and redundancy 
where 



rate \ - ^ over Z„ 



'n^3±2x_ n+l-2x \ 
v 8 ' 8 ) 
' n-3-2x n+l+2il 
v 8 ' 8 / 



if 71 = 1 (mod 
otherwise. 



Proof: Let n = Ah + 1 be a prime for some integer h and 
x 2 + Ay 2 with x = 1 (mod 4) its decomposition. Take the set 
{Co,C|} of two cyclotomic classes in F„. We compute the 
index p of this set as a DSS overZ„. By Theorem l2.5l we have 
p = min((0,2)4 + (2,2)4,(l,2)4+(3,2) 4 ).Ifn= 1 (mod 8), 
by plugging the actual values of the cyclotomic numbers of 
order four [25], we have 

(0,2)4 + (2, 2) 4 = 2(0,2) 4 



and 



S 



(1,2)4 + (3, 2) 4 = 2(1, 2) 4 



Similarly, if n = 5 (mod 

(0,2)4 + a 



and 



(1,2)4 + (3, 2) 



n + 


1 - 


2x 


have 


8 




11 — 


3- 


2x 




8 




n + 


1 + 


2x 


8 



Each cyclotomic class contains h = elements of F„. 

Hence, {Cq, C|} forms a DSS of desired parameters. ■ 
Theorem 2.9: Let n = 1 (mod 6) be a prime and binary 
quadratic form n = x 2 + 3y 2 with x = 1 (mod 3) its 
decomposition. Then the set {Cq,C|} of two cyclotomic 
classes in F„ forms a regular DSS of index p and redundancy 
rate i — 4- over Z„ , where 



-2x \ 



■ I n— 5+4x n+1- 
1111 V 18 ' 18 I 

if 2 is a cubic residue (mod n), 

■ ( n-5+ix+6y n+l-2x-12y n+l-2x+6y <\ 



18 



18 



otherwise. 



Proof: Let n = 6h + 1 be a prime for some integer h and 
x 2 + 3y 2 with x = 1 (mod 3) its decomposition. Take the set 
{Cq , Cf } of two cyclotomic classes in F„ to construct a DSS 
of index p over Z„. By Theorem 12.51 we have 



p = min((0, 3) 6 + (3, 3) 6 , (1, 3) 6 + (4, 3) a> (2, 3) 6 + (5, 3) e 



As in the proof of Theorem 12.81 a routine computation shows 
that 



(0,3) 6 + (3,3) 6 
that 

(l,3) 6 + (4,3) 6 = 
and that 

(2, 3) 6 + (5, 3) = 



18 

71 — 5+4a:+6jy 





18 


n+l- 


-2.c 


18 




n+1- 


-2x-12y 




18 


n+1- 


■2x 


18 




n+1- 


-2x+6y 



if 2 is a cubic residue, 
otherwise, 



if 2 is a cubic residue, 
otherwise, 



if 2 is a cubic residue, 
otherwise. 



Each cyclotomic class contains h = elements of F„. 

Hence, we obtain a DSS of desired parameters. ■ 
Note that Theorems 12.61 and 12.71 are special cases of these 
two classes in which the resulting DSSs are simultaneously 
perfect and regular. In general, Theorem 12.51 produces an 
optimal DSS or one close to optimal when the cyclotomic 
classes in ¥ p can be taken such that (i, e)2e + (i + e, e)^ 
are uniform or almost uniform across i. For instance, take 
n = 37 = l 2 + 4 ■ 3 2 . Then by Theorem|2j] we obtain a DSS 
of index 4 and redundancy 18 over Z37. By the Levenshtein 
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bound, the redundancy of any DSS of the same index with 
exactly two sets over Z37 must be at least 



^/2-4 - (37- 1) = 17 



which is very close to 18. 

If lower time slot usage is desirable, in principle, regular 
DSSs of better redundancy rate can be obtained in the same 
way in exchange for poorer indices p by applying cyclotomic 
numbers of higher orders. If higher indices are required 
to tolerate a higher noise level, such perfect DSSs can be 
constructed, albeit with more complicated computation, by 
taking unions of cyclotomic classes and increasing time slot 
usage accordingly (see [27 1). 

DSSs of index larger than two can be obtained by known 
recursive constructions as well. The following is a relevant 
special case of the recursive constructions given in [15]. 

— and n' 



Theorem 2.10 (hT5]l): Let n 



q-l g-1 ' 

where q is a prime power and t a positive integer. If there 
exists a DSS {Do, Di} of index p over Z„, then there exist a 
DSS {D' Q , D[} of index p 1 over Z n r, where 



\D' 



= <7 t+1 |A:| fori = 0,1 



and 



p' = mm ( W *+ 1 ) 2( 9 -l)|£>o||A|), 
and a DSS {Dg,D'{} of index p" over Z n i, where 
\D'{\ = q t+1 \D,\+in fori = 0,1 

and 

p" = min (pq t+1 ,2(q - l)|A)||£>l| + 2|D |) . 

This recursive construction gives DSSs of improved index 
by increasing the time slot usage for synchronization. For 
instance, if we apply the first half of Theorem 12.101 to the 
optimal DSS of index one over Z n with n = - q _ 1 obtained 
by Theorem 12.31 the index of the resulting DSS is at least 
qt+i _ g because we have 

2(<?- 1)7071 >2(g-l)7f 
> q t+1 - q- 

While we have focused on theoretical aspects of DSSs 
and systematic constructions, one may also look for optimal 
DSSs of specific parameters through computer searches. An 
algorithm for finding optimal DSSs is proposed in [28 1. An 
explicit example of optimal DSSs of index p over Z„ for each 
p < 5 and n < 30 can be found in ||29l . The computer 
search results and the existence of systematic constructions 
for particular parameters seem to suggest that while it is quite 
difficult to give explicit constructions, the redundancies of 
optimal DSSs of index p over Z n are generally very close 
or equal to \j2p(n — 1). 



B. Modulation layer 

We now turn our attention to the modulation layer of our 
scheme. As we have seen in the previous subsection, we can 
employ standard PPM and its variations by exploiting the 
freely available bits given by a DSS. One major benefit of 
using a DSS is that it allows for error tolerant synchronization 
while still keeping the number of time slots for synchroniza- 
tion per symbol very low. One might then wish error correctio n 
for the modulation layer at the modulation stage to eliminate 
the need of or reduce the burden on the shoulders of error 
correction at a higher level. 

Expurgated PPM is the error-correcting variant of PPM 
which can be understood as a binary constant-weight error- 
correcting code. To take advantage of our coding theoretic 
framework, we first briefly review this modulation technique 
and describe it in the language of constant- weight codes. 
A more general error-correcting variant of expurgated PPM 
will then be developed to accommodate a larger number of 
symbols. 

Expurgated PPM employs special combinatorial designs 
with cyclic automorphisms. A simple 2-design of order v, 
block size k, and index p is an ordered pair (V, B), where 

V is a finite set of cardinality v and B is a set of fc-subsets of 

V such that each pair of elements of V is included in exactly 
p element of B. Elements of V are called points while those 
of B are blocks. The simple 2-design is said to be symmetric 
if \B\ = \V\ = v. It is cyclic if the cyclic group of order v 
acts regularly on the points. 

A difference set of order v and index p is a set B of non- 
negative integers less than v such that every element of Z„\{0} 
appears exactly p times each as the difference a — b (mod v) 
between two distinct elements a, b £ B, a 7^ b. To avoid the 
triviality, we assume that k > p. Let ir be the map B M> 
B + 1 = {b + 1 (mod v) | b £ B}. It is straightforward to 
see that if the orbit Orb z „(B) = (J igZ {^(B)} is of length 
v, then the v subsets of Z v form a cyclic simple 2-design 
of order v and index p that is symmetric. Its block size is 

= 1 + ^4(^-1)^+1 

Expurgated PPM employs the v blocks of a simple 2-design 

constructed from a difference set B of order v. Trivially, the 

block set B = Orbz v (B) of this symmetric design forms 

the set of supports of the codewords of a binary constant- 

■ 1 ii-i 1 1 -it l+\/4(u— llu+1 
weight code of length v and weight k = — v 2 — 

with v codewords. Because every nonzero difference appears 

exactly p times in B, the minimum distance is 2(k — p). By 

employing this constant-weight code for modulation based on 

pulse positions as in PPM, we have Q = v tim e slots for 

each symbol interval in which K — k = 1+ V 4 ^ 
pulses are transmitted to represent M = v = Q symbols. 
Expurgated PPM thus supports the same number of symbols 
as standard PPM and has an increased minimum distance. The 
code represented by B has the property that every cyclic shift 
of a codeword is also a codeword. Hence, implementation on 
the receiver side only requires a simple correlation receiver [4| 
(see also 1*301 for implementation of signal modulation based 
on cyclic 2-designs). The above coding theoretic interpretation 
may be summarized by the following proposition: 
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Proposition 2.11: The symbols of expurgated PPM that 
transmits K pulses in each Q slot interval are equivalent 
to a binary constant-weight code of length Q, weight K, 
and minimum distance 2K ^~ if - ) with Q codewords in which 
every cyclic shift of a codeword is also a codeword. 

Proof: The set of the Q symbols of expurgated PPM is a 
subset V' K of the set V K = G | wt(u;) = if} of all 
Q-dimensional binary vectors Vi of weight K such that the set 
B = {supp(t>i) | Vi 6 T" K } of supports of the Q-dimensional 
vectors in V' K forms the block set of a cyclic simple 2-design 
of order Q, block size K, and index /i that is symmetric. 

It suffices to show that the minimum distance 2(K — n) of 
the corresponding constant- weight code is equal to 2K q > ~ ] K ^ ■ 
Because every pair of points is included in exactly /i blocks 
while each of the Q blocks contains ( 2 ) pairs, we have 

"(?)-< 



which implies that 



K(K- 1) 



- 1 



Hence, we have 



2{K- n) =2 [K- 



K(K - 1) 



Q-l 
= 2K(Q - K) 
Q-l 

as desired. ■ 
We generalize expurgated PPM by taking advantage of the 
above interpretation. Our approach is to increase the number 
of orbits over Z v while keeping the minimum distance large, 
so that the number of codewords becomes a multiple of v 
rather than exactly v. 

A (v, k, A) optical orthogonal code C C W% of length v, 
weight k, and index A is a set of u-dimensional binary vectors 
of weight k such that for any c 6 C its off-peak periodic 
Hamming autocorrelations are at most A and for any pair of 
distinct codewords c, c' e C their periodic Hamming cross- 
correlations are at most A. In other words, it is a set of v- 
dimensional vectors with k Is and v — k Os whose coordinates 
are indexed by Z v such that 



2J c t c t+ i < A 



0<t<-u-l 



for any c = (c , ci, . . . , c„_i) € C and any nonzero element 
i€Z„ and such that 



0<t<u-l 

for any pair of distinct vectors c = (co, cy, . . . , Cv-i), c' = 
(cqjC'j, . . . ,cf v _ 1 ) G C and any i G Z„. We do not consider 
the trivial case k = A and always assume that k > A. 

Optical orthogonal codes have been extensively investigated 
from various viewpoints including the initial motivation in the 
context of code-division multiple-access fiber optical commu- 
nications [31]. A useful observation for our goal is that an 
optical orthogonal code is equivalent to a constant-weight code 



in which every cyclic shift of a codeword is also a codeword, 
which is the property we would like for our signal modulation 
purpose. To see the equivalence, take the union T> of a (v, k, A) 
optical orthogonal code C and the set of all cyclic shifts of each 
codeword. Then T) forms a constant-weight code of length v, 
weight k, and minimum distance at least 2(k — A). Trivially, 
the converse also holds. Because we included all cyclic shifts 
of c € C in T>, every cyclic shift of d e D is naturally in 
T> again. Because the condition that k > A is assumed, the 
number of codewords of the corresponding constant-weight 
code is v\C\. 

Now by Proposition 12.111 if we look at expurgated PPM 
with v symbols that uses k pulses per symbol in our coding 
theoretic framework, it is a constant-weight code of length 
v, weight k, and minimum distance 2(k — A) with exactly v 
codewords in which every cyclic shift of a codeword is also a 
codeword. In other words, it is simply the union of a (v, k, A) 
optical orthogonal code with only one codeword and its v — 1 
cyclic shifts. Because it is also a symmetric design, we have 
A = Mtil. Thus, we obtain the following proposition: 

Proposition 2.12: The set of symbols of expurgated PPM 
that transmits K pulses in each Q-slot interval are equivalent 
to a (Q, K, K q~-P ) optical orthogonal code with exactly one 
codeword. 

Since an optical orthogonal code C of length v gives rise to 
a constant- weight code with v\C\ codewords by joining the v 
cyclic shifts of all c G C, it is natural to generalize the PPM 
method such that modulation exploits an optical orthogonal 
code with more than one codeword. With this generalization, 
a larger number of symbols can be supported compared to 
expurgated PPM while still using the same decoder for each 
orbit and maintaining the error correction mechanism at the 
modulation stage. Summarizing the discussion given above in 
this subsection, we have the following theorem: 

Theorem 2.13: The set of the codewords of a (v, k, A) op- 
tical orthogonal code C and all cyclic shifts of each codeword 
defines a (u|C|)-ary signal modulation technique in which each 
of the v\C\ symbols with mutual Hamming distance at least 
2{k — A) is represented by k single pulses transmitted at k out 
of v time slots. 

We would like optical orthogonal codes with the largest 
possible number of codewords for given v, k, and A. Because 
an optical orthogonal code C is already a constant-weight code 
before joining the cyclic shifts of each codeword, the Johnson 
bound gives the upper bound on the number of codewords of 
C: 

Theorem 2.14 filJ2l/): Let C be a (v, k, A) optical orthogo- 
nal code. Then it holds that 

v- A 



\C\< 



v-1 
k-l 



v-2 



k-X 



A (v,k,X) optical orthogonal code is optimal if the number 
of codewords attains this upper bound. Table [II] lists well- 
known classes of optimal optical orthogonal codes of index 
one with more than one codeword that span a variety of 
lengths and weights. There are numerous other constructions 
and existence results. For the case when the index is one, all 
other known results can be found in ll42~l -[47 l and references 
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TABLE II 

Some classes of optimal optical orthogonal codes of index one with more than one codeword 



I Ulllltll ft 1 


VVclglll ft 


1 11 f 1 ov \ 

lllUtA A 




| 1 Aflcfl , 'lllli 


tvLIC I cIlCc 


n 


3 


1 




n ^ 14, zU (inocl z4 j 


11 1 1 1 111 


n 


1 


1 




n = 0,6,18 (mod 24) 
p = 1 (mod fc(fc — 1)) prime, 


03) 03 


V 


any 


1 


L*(*-i)J 


V > Ck, 
Cfc constant dependent on k b 


03 




q 


1 


q-1 


q prime power 


Affine geometry with origine deleted 


9-1 


9 + 1 


1 


— ri I even 
Ife *°dd 


q prime power 


Projective geometry 



d This is a necessary and sufficient condition for the existence (see 1381 for a short proof). The constrains for the other classes are sufficient conditions. 

b C4 — C5 — 1391 . cq — 61 1401 . For k > 7 in general, — 1 '- Note that if p < 2k(k — 1) + 1, the code does not have more than one codeword. 

L The same parameters may be realized as a generalized Bose-Chowla family 1411 . 



therein. For the latest results on optical orthogonal codes 
of higher index, we refer the reader to [|42l . [4-8 — [50 1 and 
references therein. 

III. Concluding remarks 

We introduced a coding theoretic framework to the study of 
signal modulation based on pulse positions in the time domain. 
With the further help of combinatorial design theory, this 
approach allowed us to develop a self-synchronizing scheme 
with significantly improved efficiency while maintaining com- 
patibility to existence modulation techniques such as standard 
pulse position modulation. In fact, the number of time slots 
required per symbol for synchronization is now improved 
from f2(Q) to 0(Qz) in the asymptotic sense, breaking the 
fundamental limit of the previous synchronization method. 

We were also able to generalize the recently introduced 
error-correcting pulse position modulation technique to realize 
a larger number of symbols. This generalization is particularly 
appealing as the modulation layer for our synchronization 
method because this way error correction can be fully sup- 
ported at the modulation stage. 

In the previous section on synchronization and error cor- 
rection, we have placed particular emphasis on generality so 
as not to unnecessarily spoil the potential. For this reason, 
our focus has been on the minimum distance. However, it 
would also be of importance to investigate various other 
aspects of our scheme by assuming a particular context. In 
fact, PPM and its variants have extensively been studied with 
various applications in mind (see, for example, J5T), J52 | for 
use in ultra-wideband communications, ll53l for free-space 
optics scenarios, and [54] for the purpose of atmospheric 
optical code-division multiple-access). Among many possible 
directions of more focused research, it would be of particular 
interest to analyze the spectral efficiency in detail and more 
accurately estimate the synchronization error rate and bit error 
rate over a reasonably realistic channel. 

If one wishes to extract finer structural information than 
minimum distance for a detailed analysis in a specific context, 
the algebraic properties of DSSs and optical orthogonal codes 



may be effectively exploited. For instance, it is straightforward 
to see that if we employ an optical orthogonal code C of length 
v, block size k, and index 1 to form the constant-weight code 
T> for modulation by joining cyclic shifts, for any codeword 
d e T> the number nd of codewords of Hamming distance 
2 (k — 1) from d, which are the nearest, and the number fa of 
other codewords except d, which are all of Hamming distance 
2k from d, are 

rid = k \C\ — k 

and 

f d = v\C\ - 1 -rid 
= {v - k 2 )\C\ + k - 1 

respectively. 

Another aspect we did not address is the relation of our 
version of expurgated PPM to error-correcting codes at a 
higher level. If we see PPM as a constant-weight code, the idea 
of error correction at a higher level is a code concatenation 
in one sense. There have been proposed various types of code 
and use for such concatenations (see, for instance, |55|-[59|). 
It would be of interest and importance to understand how best 
to exploit our scheme along with an error-correcting code at 
a higher level in a particular communications system. 
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